Effect of realistic finite-size impurities on T c in Fe-based superconductors 
based on the five-orbital tight-binding model 
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We present a systematic study of the impurity effect on T c in Fe-based superconductors, assuming 
that the sign-reversal s-wave state due to inter-pocket repulsion (s±-wave state) is realized. For this 
purpose, we introduce several realistic impurity models with non-local modifications of potentials 
and hopping integrals around the impurity site. When we use the impurity model parameters for 
3d- and 4d-impurity atoms derived from the recent first principle study by Nakamura et at, we find 
that the s±-wave state is very fragile against impurities: The superconductivity without impurities 
Tco = 30K is destroyed by introducing small residual resistivity p" = 5z _1 ~ lO-z^ 1 [/iQcm] 
(z^ 1 = m*/m being the mass-enhancement factor), consistently with the previous theoretical study 
for the on-site impurity model by Onari and Kontani. This result is essentially unchanged for 
different non-local impurity models with realistic parameters. We also discuss the effect of the 
impurity-induced non-local orbital order on the superconducting state. 

PACS numbers: 74.20.-z, 74.20.Fg, 74.20.Rp 



I. INTRODUCTION 

Since the discovery of Fe-based high-T c superconduc- 
tors [1] , the symmetry and the gap structure of the super- 
conducting (SC) state have been studied very intensively. 
It had been established experimentally that s-wave (A\ g 
symmetry) SC state is realized in usual Fe-based super- 
conductors. The gap structure in many optimally-doped 
high-T c compounds is nearly isotropic and fully-gapped 
[2-4], although some compounds show accidental nodal 
gap structure. In the phase diagram, the SC phase is 
realized next to the orthorhombic structure transition at 
Ts, and the magnetic order is also realized at Tn < T$. 
Below Ts, the orbital polarization n xz ^ n yz is realized 
[5], and sizable softening of shear modulus Cqq [6-8] in- 
dicates the development of orbital fluctuations near the 
orthorhombic phase. Strong spin fluctuations are also 
observed near the magnetic ordered phase. 

These observed orbital and spin quantum criticalities 
have been intensively studied theoretically, since they 
would be closely connected to the pairing mechanism. 
Within the random-phase-approximation (RPA), strong 
spin fluctuations develop in the multiorbital Hubbard 
models for Fe-based superconductors. Therefore, spin- 
fluctuation-mediated s-wave state with sign reversal (s±- 
wave state) is obtained by the RPA [9-12]. However, the 
RPA fails to explain the non-magnetic structure tran- 
sition at Ts. Also, orbital-fluctuation-mediated s-wave 
state without sign reversal (s ++ -wave state) is realized 
by introducing the quadrupole interaction <? qua d due to 
Fe-ion oscillations [13-15]. Even for <? qua d = 0, strong or- 
bital fluctuations are obtained by improving the RPA by 
including the vertex correction (VC) for the susceptibility 
that is dropped in the RPA [16, 17]: Since spin and or- 
bital fluctuations mutually develop in the self-consistent 
VC method, both s±-wave and s++-wave states can be 
obtained by solving the multiorbital Hubbard model. 



To distinguish between these two SC gap states, vari- 
ous phase-sensitive experiments had been performed [18- 
25]. For example, inelastic neutron scattering experi- 
ments had been performed to find the magnetic resonance 
scattering due to the sign reversal [19-21]. However, ob- 
served "resonance-like" hump structure can be explained 
even if s++-wave state is realized [26]. Impurity effect 
measurement is another significant phase-sensitive exper- 
iment since T c would be strongly suppressed by inter- 
pocket impurity scattering of Cooper pairs if s±-wave 
state is realized [22-25]. Many theoretical studies have 
been devoted so far [12, 27-32]. In Ref. [30], the present 
authors studied the impurity effect due to local nonmag- 
netic impurities based on the realistic multiorbital model, 
and showed that T c in the s±-wave state is strongly sup- 
pressed by inter-pocket impurity scattering of Cooper 
pairs. However, effect of the possible non-locality of the 
impurity potential had not been studied. 

Therefore, in this paper, we present a quantitative 
study of the non-local impurity effect on the s±-wave 
state based on the realistic five-orbital model for Fe- 
based superconductors. We introduce realistic models 
of Fe-site substitutional impurity atoms, which contains 
the non-local modifications of potentials {I,V\,V2) and 
hopping integrals (At) around the impurity site, refer- 
ring to a recent first principle study by Nakamura et 
al. [33]. In various non-local impurity models with 
realistic parameters, the s±-wave state is very fragile 
against impurities: Original transition temperature at 
T c0 = 30K is destroyed when the residual resistivity is 
just pjf = 5z _1 r~j lOz -1 [/xf2cm], where = m*/m is 
the mass-enhancement factor. This result is consistent 
with the previous theoretical study for the on-site impu- 
rity model in Ref. [30]. Thus, experimentally observed 
robustness of T c against impurities in various (optimally- 
doped) Fe-based superconductors [22-25] would indicate 
the realization of the s+_|_-wave state. 
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The effect of Fe-site substitution on T c had been care- 
fully studied in La(Fe,A)As(0,F) for A=Co,Ni [34] and 
in Ba(Fe,X)2As2 as well as Ba(Fei_ x _ s/ Co a: Cuy)2As2 
[35-37]. In these materials, T c is well scaled by the 
amount of carrier doping caused by A-atom substitu- 
tions, not by the impurity concentration irrespective of 
the large impurity potential [37]. This fact would mean 
the robustness of the SC state against strong impurity 
scattering in these materials as stressed in Ref. [34]. 

In Fe-based superconductors, impurity potential ma- 
trix I\,(k,k') in the band-diagonal basis is fc-dependent, 
reflecting the multiorbital band structure [30]. How- 
ever, their fc-dependences had been frequently neglected 
in previous studies for simplicity. In this "constant Jb 
model" , both the intra-band and inter-band scatterings, 
Jb and J b , are constant parameters. In this model, the 
s_l--wave state can be stable against impurities when 
l^b/^b| "C 1, although it seems unrealistic since both 
hole-Fermi surfaces (h-FSs) and electron-Fermi surfaces 
(e-FSs) are composed of the common c£-orbitals [30]. In 
addition, obtained results are strongly changed in the 
unitary and intermediate regimes once the fc-dependence 
of Jb is taken into account; see Sec. VC. 



= WrW&aCfV*, (2) 

r,r' 1,1' ,a 

where J, V\ and V2 are the on-site, nearest-neighbor 
(NN), and next-nearest-neighbor (NNN) impurity poten- 
tial, and &toirl> is the modulation of the NN or NNN 
hopping integrals between site and site r. The present 
impurity potential model is depicted in Fig. 1. 
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FIG. 1: (Color online) The present non-local impurity poten- 
tial. 



II. MODEL HAMILTONIAN 

In this paper, we study the impurity effect based on the 
realistic two-dimensional five-orbital tight-binding model 
[9]: 



H 



k.a.l.l' 



(i) 



where h k is the 5x5 matrix given by the Fourier trans- 
formation of the hopping integral r , v introduced in 
Ref. [9]. Here, 1,1' represents the orbital indices, and a 
is the spin index. The matrix elements of h k is given by 
the Fourier transformation of the hopping integral r , v . 
When the electron filling per Fe-site is n — 6.0, There 
are two hole-pockets around the T point, one hole-pocket 
around (71", 71") point, and two electron-pockets around 
(n, 0) and (0, tt) points. 

In addition, we introduce the following nonmagnetic 
and non-local impurity potential at site = (0, 0): 
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III. GAP EQUATION, CALCULATION OF 
RESIDUAL RESISTIVITY 

In the present model, the T-matrix due to infinite num- 
ber of impurity scattering processes is given as 



T[ie n ) = W(l — g(ie n ) ■ W)~ 



(3) 



where g r iyv is the free Green function in real space, 
given by the Fourier transformation of G k (ie n ) = (ie n + 
M — an d e n = (2n+l)7rT is the fermion Matsubara 

frequency. 

When the impurity concentration is dilute (n- uup C 1), 
the normal self-energy above T c is well approximated by 
the T-matrix approximation. It is given by 



5T, k (ie n ) = n imp f kjk ,(ie n ) 



k=k' 



(4) 



where T k j~>(ie n ) is given by the Fourier transformation 
of eq. (3). The retarded (advanced) self-energy is 
given by the analytic continuation ie n —r e + iS ■ sgn(n): 

5±l R{A \e) = St,(k,e + (-)iS). Then, the Green function 
in the band-diagonal basis is G^ k (e) = l/(e + fj, — E ak — 
27afc(e)); where a is the band basis, E ak is the dispersion 
of the ath band, and 7 Q fc(e) = — Im<5S£ R (e) is the quasi- 
particle damping rate due to impurities. 7 a fe(0) can be 
rewritten as 



''imp 
47T 



X 



-? JFS/3 V(3k> 



IT&^OF, (5) 
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where J FSa dk is the integration on the FSa. 

When the impurity concentration is low enough, the 
inter-band contribution to the conductivity is negligible. 
If we drop the current vertex correction (CVC), the con- 
ductivity is given by 



band 



no CVC 



k , a 



de 



de 



\Gl k {e)\\l, v , (6) 



where v = x or y, /(e) = (e £/T + 1) -1 , v a k,v = dE ak /dk v 
is the quasiparticle velocity, and c is the inter-layer dis- 
tance. However, to obtain the exact conductivity for 
«imp <C 1, the CVC should be taken into account. The 
exact expression for the conductivity is given as 

<^ = 7 E / 7 (-%) IG^WI 2 w JohAcW) 

k,a J \ ' E ah 

Hereafter, we put c = 0.6nm that corresponds to Bal22 
compounds. J a k. v is the total velocity with CVC, which 
is given by solving the following Bethe-Salpeter equation: 

+n imp J2 \ T *k,W (e)| 2 |G£ fe ,(e)| 2 J/3 fe >(e), (8) 
k',0 

where T^ k /3fc , (e) is the T-matrix in the band-diagonal 
basis. Here, we neglect the inter-band contribution since 
it is negligible when n; mp <C 1. 

Next, we study the impurity effect on T c . The gap 
equation at T = T c in the band basis is given by 

Z ak {^n)A ak {e n ) = __£.V^ f V a k ,0fc'(e n , e m ) 

4n 7~i JFS/3 v l3k> 



X^^+SKk^n), 



(9) 



where A ak is the SC gap function in the band-diagonal 
basis, and V is the pairing interaction. Z ak (e n ) = 1 + 
7afc(e™)/l e n|, and <5£^ is the linearized anomalous self- 
energy is given as 

5Y> a ak {ie n ) = n imp ^\T ak ^ k ,(ie n )\ 2 \Gpk'(i^n)\ 2 ^f)k'(^n 

fe',/3 

dk ' ]rr (■ M2 A /3fc'(» £ ") 
\J-ak,0k'^£n)\ j j : 



E 



4?r y JFS/3 u /3fe' 



T c against impurities is derived from this relation, known 
as the Anderson theorem. However, Anderson theorem 
is totally violated in the s±-wave state. 

In the present study, we consider the case of spin fluc- 
tuation mediated intra-orbital repulsive interaction be- 
tween e-FS and h-FS. Then, the spin susceptibility is ap- 
proximately diagonal with respect to the d-orbital, since 
it is strongly enhanced by the intra-orbital repulsion U. 
That is, spin-spin correlation develops only in the same 
d orbital. By taking this fact into account, we introduce 
the following repulsive pairing interaction: 



(10) 



By solving the gap equation (9) together with eqs. (4) 
and (10), the impurity effect on T c for n; mp <C 1 is exactly 
calculated within the BCS theory. 

The reduction in T c is caused by the function ~f ak ■ 
Aafc/knl — <5S° fc , which vanishes if A ak is independent 
of band and momentum (^isotropic s++-wave state) ac- 
cording to eqs. (5) and (10). Then, the independence of 



d— orbital 



V a k,{Sk<(e n ,e n <) = 5 E \ u iA k )\ 2 \ u iA k ') 
i 

xD(e n )D(e n »), 



(11) 



for a <Ge-FS and (3 eh-FS and vise versa, and V a k,j3k f = 
for a, /3 ee-FS or a,/3 eh-FS. Ui >a {k) = (l,k\a,k) is 
the unitary matrix connecting between band-basis and 
orbital-basis, and D(x) = w^/(x 2 +w 2 ), where uj c is the 
cut-off energy of the pairing interaction. In the present 
study, we put g = 3.22eV and w c = 0.03eV. In this case, 
a fully-gapped s±-wave state with T c o = 30K is realized 
when rii 

eV, and the unit of temperature is K. 



imp = 0. Hereafter, we set the unit of energy is 



IV. NUMERICAL RESULTS 
A. Nakamura's impurity model 

First, we study the non-local impurity potential 
W rl,r'l' for A-atom (A= Mn, Co, Ni, Zn, Ru) given by 
the first principle study by Nakamura et al. [33] . The de- 
rived values of the on-site impurity potential / are +0.28, 
-0.35, -0.87, -8.05 and -0.02 (the unit is eV) for Mn, 
Co, Ni, Zn and Ru, respectively, as shown in Table I. 
Thus, Zn atom works as the unitary impurity scattering 
center, whereas other 3d atoms (Mn, Co and Ni) induce 
intermediate impurity scattering. The off-site impurity 
potentials (Vi and V2) are very small. In addition, strong 
modulations of NN and NNN hopping integrals (A^ 1 ) 
and At' 2 ') are induced around the impurity atom. Espe- 
cially, At^ /t° for Ad atom Ru takes large positive value, 
reflecting the larger radius of 4<i-orbitals. Thanks to the 
Nakamura's model, we can present a quantitative analy- 
sis of the impurity effect on the s±-wave state, by taking 
the non-local nature of the impurity potential appropri- 
ately. 

Figure 2 show the obtained T c as function of the (a) 
impurity concentration ni mp and (b) residual resistivity 
po = 1/<7 X for various impurity atoms in the case of 
= 1. Although the Fe-site substitution induces the 
"impurity potential" and "carrier doping" , we neglect the 
latter effect by fixing the electron filling n = 6.0, in order 
to concentrate on the former effect. When the mass- 
enhancement factors is finite, the reduction of T c 
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Mn 


Co 


Ni 


Zn 


Ru 


I Ave 


+0.28 


-0.35 


-0.87 


-8.05 


-0.02 


A/V(/ A vc) 


-1.14 


+1.09 


+2.05 


+3.91 


-0.02 



TABLE I: 7a vc is the averaged on-site impurity potential in 
Nakamura's model [33]. The unit of energy is eV. AN (I Ave) 
is the change in the local electron density at the impurity 
site due to Tavg in the present model without interaction. In 
Nakamura's model, the on-site potential is orbital-dependent, 
and At w and At (2) are included. 



Po( n imp) are multiplied 
(a), the critical impurity 



W x } 



per impurity concentration, — (T c — T c o)/ni mp , is renor- 
malized by z, while po/^imp is independent of z [38]. 
Therefore, both n? mp and p" 
by z _1 . According to Fig. 2 
concentration for the disappearance of T c , nf^ pl strongly 
depends on the impurity atoms: nf* = 0.6z _1 [%] ~ 
3.5z _1 [%] for 3d-impurities (Mn, Co, Ni, Zn), while 
nf mp = 24z _1 [%] for Ru-impurities. In contrast, the 
values of p" shown in Fig. 2 (b) are almost independent 
of impurities for 3d- impurity atoms (~ 5z _1 [//Ocm.]), 
while p" ~ 10z _1 [/ificm] for Ru-impurity. 

The residual resistivity given by Nakamura parameter 
is very small (po/^imp < 5 [pflcm/%]) expect for 
X = Zn. One of the reasons would be that the impurity 
potential given in Ref. [33] may be normalized, although 
the bare impurity potential is required for the present 
study. (The normalization would caused by the modifi- 
cation of the wavefunction around the impurity site in 
solving the Kohn-Sham equation, while this process is 
also included in the T-matrix.) By taking this effect into 
account, we study the impurity potential M x W^ lr , v 
with M > 1. To see the effect of M, we analyze the 
case of M = 4 in Fig. 3 for instance. We show (a) 
"■imp-dependence and (b) po-dependence of T c . Note that 
similar results are obtained for M > 4. In this case, 
T c q = 30K is suppressed by a few percent impurity con- 
centration. In contrast, p§ r is approximately independent 
of M for all impurity atoms. We stress that po /ni mp is en- 
larged to ~ 50 [pQcm/%\ by taking the impurity-induced 
non-local orbital order into account, because of the en- 
larged impurity scattering cross section [39] . In this case, 
s±-wave state is strongly suppressed, as we will show in 
Sec. IV C. 

For the convenience of analysis, we hereafter study the 
ratio R = -(T c -T c0 )/po = -AT c /p at n imp « < lp : R 
is rather independent of the impurity potential strength, 
and also R is essentially independent of T c q accord- 
ing to the Abrikosov-Gorkov theory. In Fig. 4, we 
show the obtained R as function of M for various im- 
purity atoms. As recognized in Ref. 2, the relation 
Too/p" ~ 1.5 Ji is satisfied for all atoms. By taking the 
mass-renormalization factor into account, we obtain that 
R = 3z ~ 5z [K///f2cru] for all 3c?-impurity atoms studied 
in Ref. [33], except for the Ru-impurity. (The horizon- 
tal broken line Ri =oa — 3.6 [K///f2cm] is the value for 
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FIG. 2: (Color online) Obtained T c as function of (a) m m p 
and (b) po using the impurity potential W£[r'l' f° r ^-atom 
(X= Mn, Co, Ni, Zn, Ru). 



the infinite on-site impurity potential (/ = oo) studied in 
Ref. [30].) 

In contrast, experimentally observed R in optimally- 
doped 1111 compounds (z _1 ~ 2) and 122 compounds 
(z' 1 - 3) is i? cxp - 0.1 [K/pOcm] [22-25]. Therefore, 
the s±-wave state would be too fragile against nonmag- 
netic impurities to explain experimental robustness of T c 
against impurities. Similar result was reported by Ikeda 
et al [40] by using the Nakamura's impurity model. 



B. ideal non-local impurity models 

According to the first principle studies in Refs. [33, 41, 
42] , the impurity-induced change in the electron density 
is strongly localized at the impurity center. This fact in- 
dicates the smallness of the non-local impurity potentials 
Vi in eq. (2). However, to obtain useful knowledge on the 
impurity effect, we introduce four ideal impurity poten- 
tial models shown in Fig. 5 (a): (i) /-model (only on-site 
potential), (ii) V\ -model (only NN potential, without J), 
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FIG. 3: (Color online) Obtained T c as function of (a) riimp 
and (b) po using the impurity potential 4 x r , v . 

(iii) V2-model (only NNN potential, without /), and (iv) 
^As-model (plaquette impurity potential due to As-site 
substitution). Note that (ii) and (iii) are very unrealistic 
potentials. Here, we put Av-*' — for simplicity. 

Figure 5 (b) shows the ratio R = — (T c — T c0 )/p at 
Himp = 0.1 [%] for models (i)-(iv) as function of the impu- 
rity potential, in the case of z^ 1 — 1. The obtained R for 
the V2-model is as large as that for /-model, while those 
for other two models are smaller when Vi, Va s < 0. Espe- 
cially, R < 1.0 [K//if2cm] is realized for the Vi-model for 
— 1 < Vi < 0.5. We also show —AT c /rii mp and po/^imp 
with CVC in Figs. 5 (c) and (d), respectively. Com- 
pared to the /-model, — AT c /n; mp in the V± model is 
comparable, while p in the V± model is much larger. 
For this reason, the relation model ^ Rvi— model can 
be achieved. 

However, the above Vi-model without on-site potential 
/ is very unrealistic. Thus, we also introduce / to the 
Vi-model: Figure 6 shows the obtained (a) R and (b) 
Po/^imp as function of V± with finite /. In the case of 
|/| > 0.5, —AT c /p" quickly approaches to the value for 
the /-model for | V\ | < 0.2. In real compound, the relation 
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FIG. 4: (Color online) R = -(T c - T c0 )/po at n imp = 0.1 [%] 
for the impurity potential M x W x for various X atoms. The 
horizontal broken like at Ri=oa = 3.6 [K/pf2cm] represents 
the value for the infinite on-site impurity potential (/ = oo) 
studied in Ref. [30] for z~ x = 1. 

| Vi | < 0.1 is expected, as we will discuss in Sec. VA. 

We also study the effect of hopping integral inhomo- 
geneity around the impurity site. For a systematic study, 
we define At^ in eq. (2) as Atgi rl , = x\ ■ t u; rV for the 
NN sites (\r\ = 1), and At { ^ rl , = x 2 -t° l rll , for the NNN 

sites (|r'| = \2). In the case of Ru-impurity, X\ and x 2 
are positive according to Ref. [33] , because of the larger 
radius of 4<i-orbitals. However, we also study the case 
Xi < since this situation might be realized by irradia- 
tions, by shifting the Fe-ion position outside of the FeAs 
plane. 

In Fig. 7, we show the obtained R and po/^imp for 
various on-site potential / as function of x%, in the case 
of (a)(b) X2 = and (c)(d) X2 — x\. In the former 
case, the obtained R is very small for / = 0, while it 
quickly increases for finite /. In contrast, R in the latter 
case is large even for / = 0. In both cases, the residual 
resistivity is very small for x\ > 0, since the magnitude 
of the hopping integral is increased locally. Therefore, 
the s±-wave state is strongly suppressed by the hopping 
integral inhomogeneity in both xi = and X2 = x\ cases, 
except for a special case / = x% = 0. 



C. impurity-induced non-local orbital order 

To study the effect of a very wide-range impurity po- 
tential, we analyze the effect of the impurity-induced 
non-local orbital order (NL-OO) derived in Ref. [39] 
on the s±-wave state. When this impurity-induced NL- 
OO is formed, the residual resistivity per 1% impurity 
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FIG. 5: (Color online) (a) Impurity potentials for the /-model, 
Vi-model (without /), V2-model (without /), and VAs-model. 
(b) Obtained R = — (T c — T c o)/po as function of the impurity 
potential. Small R for Vi-model is enlarged by introducing the 
on-site potential / if |/| 3> |Vi|; see Fig. 6. (c) — AT c /m m p at 
riimp =0.1 [%] and (d) po/?iimp as function of the impurity 
potential. 



atoms increases to ~ 50 |/tOcm] , because of the enlarged 
impurity scattering cross section [39]. This fact would 
resolve the problem that the residual resistivity derived 
from Nakamura parameter is very small expect for Zn- 
impurity atom. 

Recent discovery of "electronic nematic transition" in 
the tetragonal phase, free from any lattice deformation, 
has been attracting great attention. For example, in 
"detwinned" Ba(Fei_ a; Co a ;)2As2 [43] under very small 
uniaxial pressure (~ 5MPa), sizable in-plane anisotropy 
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FIG. 6: (Color online) (a) R for the Vi-model with on-site 
potential /. For finite |/| (> 0.5), R ~ i?/ =00 = 3.6 [K//if2cm] 
is realized for |Vi| ~ 0.1. (b) The residual resistivity po/^imp. 
The value Po = °°/ n im P = 10.6 [/ificm/%] is shown by broken 
line. 




FIG. 7: (Color online) Obtained R and po/nimp as function of 
xi — A^ 1 ' /t°, in the case of (a)(b) = and (c)(d) xi = X\. 
R ~ Ri=oo = 3.6 [K/fiUcm] is realized for |/| > 0.5. 



of resistivity emerges at T*, which is about 10K^100K 
higher than T§. The nematic order is also observed in 
BaFe2(As,P)2 by the magnetic torque measurement [44]. 
In Ref. [39], the authors discussed the impurity-induced 
electronic nematic phase using the mean-field approxi- 
mation in real space, by introducing the quadrupole in- 
teraction g qua d- When orbital fluctuations develop, local 
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impurity potential induces NL-00 with C2 symmetry, 
actually reported by STM/STS autocorrelation analyses 
[45, 46] . The large cross section of the NL-00 gives giant 
residual resistivity, When C2 nanostructures are aligned 
along a-axis, the in-plane anisotropy of resistivity reaches 
40%, consistently with experiments [43]. 



D. derivation of 171,, 
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FIG. 8: (Color online) Obtained T c of the s±-wave state in 
the presence of impurity-induced NL-OO obtained for I = — 2 
given in Ref. [39]. 

In Fig. 8, we show the obtained T c as function of 
"-imp for the (i) on-site impurity potential I = —2 and 
(ii) impurity-induced C2 orbital order obtained by the 
mean- field approximation in real space for I = —2 and 
ffquad = 0.218, shown in Fig. 1 (c) of Ref. [39]. In the case 
(ii), the critical impurity concentration is just n?^ p = 
0.3z _1 [%], which is about one-third of nj^ for the case 
(i). The corresponding in-plane averaged resistivity is 
pjf = 17.3z _1 [/ir2cm], which is about four times pg r for 



the case (i). That is, p" is enlarged while n" np is reduced 
when the impurity-induced NL-00 is realized. Note that 
the enlarged p" is still one order of magnitude smaller 
than experimental values (300 ~ 500 [//f2cm]). 

Therefore, impurity-induced C2 orbital order should 
cause strong suppression of the s±-wave SC state. Simi- 
lar behaviors (larger p" and smaller nf*) are also real- 
ized by impurity-induced short-range AF order in nearly 
AF metals, such as under-doped cuprates [47]. We con- 
sider that this impurity-induced NL-00 will also sup- 
press the s++-wave state moderately, because of the sup- 
pression in the density-of-states (pseudo-gap formation) 
due to the orbital (short-range) order. Under T c , the 
impurity-induced NL-00 will cause the "Swiss cheese 
hole state" in the s++-wave state. 



We have introduced various non-local impurity poten- 
tials, and studied the impurity effect on the s± wave state 
driven by the pairing interaction in eq. (11). We find 
the relation R = 2.5z ~ 5z [K/pflcm] holds for 3d- and 
4d-impurity atoms, whereas R — lz ~ 2z [K//if2cm] in a 
special model. To understand these numerical results, we 
analyze the averaged T-matrix between FSa and FS/3: 



IT;, 



e-FSh-FS 

:E E 



/fsq dk J FSj3 dk'\T ak 



:ik' 



-(12) 



IT:, 



f FSa dk f FSj3 dk'\T a k.j)k 12 



slE + E 



a,fi a,/3 



(13) 



Then, the averaged ratio between inter-pocket and 
intra-pocket scattering amplitude would be t 2 = 

|Ti n ter| 2 /|Ti n tra| 2 - 

Figure 9 (a) shows the values of R in various im- 
purity models as function of 1 = 2t 2 / (1 + t 2 ) = 
2|T intor | 2 /(|T intor | 2 + |T intra | 2 ). The CVC for the conduc- 
tivity is taken into account correctly. It is found that 
R is approximately scaled by x for various kinds of im- 
purities. This result is naturally understood since the 
s±-wave state is suppressed by the inter-FS scattering, 
whereas both inter- and intra-FS scattering contribute 
to po oc 7, as understood by eqs. (5) and (10). There- 
fore, the following relationships would be realize [12]: 



-AT C cx |T ntor | , 

I 2 + IT, 



PO OC |T ir _ 



(14) 
(15) 



We will discuss this issue in more detail in Sec. VB. 

Another important finding in Fig. 9 (a) is that the 
value of t = |Ti n t cr |/|T n tcr| in the five-orbital model 
is approximately independent of the impurity potential 
strength (in both Born and unitary regimes), for all the 
impurity models studied here. This fact means that the 
"constant /b-niodel can be applicable only for Born im- 
purities, as we will discuss in Sec. VC. 

In Fig. 9 (a), the values of R for V\- and VAs-niodels 
apparently deviate from the y = 3.5a: line. This fact 
originates from the large contribution from the CVC, 
which becomes important when the impurity potential 
has finite scattering cross section, since the forward im- 
purity scattering is correctly subtracted by taking the 
CVC into account. In fact, as shown in Fig. 9 (b), the 
relation R cx x is apparently improved by neglecting the 
CVC. Therefore, we should take the CVC into account 
to obtain quantitatively reliable value of R, especially for 
wide-range impurity potentials. 



V. DISCUSSIONS 

In the previous section, we calculated the impurity ef- 
fect on the s±-wave state for various non-local impurity 
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FIG. 9: (Color online) Obtained R's as function of x 



2\Tw 



7(|r in 



+ |T in 



in the case of (a) with CVC 



and (b) without CVC, for 3d- and 4<i-impurity atoms as well 
as /-model, VAs-model, and Vi-model without I: Note that 
the last potential is unrealistic. |Tinter| = |Tintra| corresponds 
to x — 1. The proportional relation R oc x becomes worse by 
including the CVC, especially for the Vi-model without / and 
VAs-model. 



larger than V\ in magnitude in real impurities. 

By introducing V\ the electron number at the NN site 
of the impurity center is changed by ATVnn ~ — 2 x V\ 
for |Vi| <C 1, according to the analysis of the present 
five-orbital tight-binding model. However, according to 
the first principle study [41], |AtVnn| would be at most 
0.1 ~ 0.2, meaning that |V1| < 0.1 due to the strong 
screening effect in real compounds. Therefore, the effect 
of V\ would be negligible in the study of the impurity 
effect in Fe-based superconductors. 



B. Why |Tinter| < |T in tra] in the Vi-model? 

We have shown in Fig. 9 that the relation |Ti nter | <C 
|21ntra| holds in some non-local impurity potentials. To 
understand the reason, we consider a orbital-less square 
lattice model for simplicity. In the Born approximation, 
the T-matrix for the Vi-model is given as 



T k<p (e) = 4^3(0,0) (e) 

+4V 1 2 (/( 1;1 )(e)(cos k x cospy + cos k y cosp x ) 
+2^3(2,0) (e)(cos(k x -p x ) +cos(k y -p y )), 



where g( x . y )(e) is the Green function in real space at 
r — (x,y). The cosine terms in eq. (16) originate from 
the non-locality of the impurity potential. For the in- 
tra hole-pocket scattering (fe = p = 0), eq. (16) gives 
Tintra = 4V 1 2 (5( 0: o) + 2 ff(i,i) + 5(2,o))- On the other hand, 
for the scattering between hole- and electron-pockets 
(k = (0,0) and p = (tt,0)), eq. (16) gives T intcr = 
4V 1 2 <?(o,o)- Therefore, the relation |Ti nter | <C |?l ntra | would 
be possible when <7 (0)0 ) - g {lX) ~ 3(2,0) ■ 



C. Comparison with the constant Jb model 



models. When we use the impurity model parameters ob- 
tained by a recent first principle study by Nakamura et 
al, the s±-wave state is fragile against impurities, consis- 
tently with the previous theoretical study for the on-site 
impurity model by Onari and Kontani [30]. This result 
is qualitatively unchanged for different non-local impu- 
rity models, when the parameters are realistic. Here, we 
present more detailed discussions. 



A. Estimation of the value of Vi 

In Fig. 5 (a), we have shown that R — — AT c /po for the 
Vi-impurity model is as small as lz [K//if2cm], which is 
still much larger than experimentally observed relations 
R < 0.1 [K/^Ocm] in single crystals [23, 24]. Moreover, 
as shown in Fig. 6 (a), R is strongly enlarged by in- 
troducing the no-site potential /, which should be much 



In this paper, we analyzed the impurity effect based 
on the realistic five-orbital model with non-local impu- 
rity potentials. On the other hand, more simple two- 
band model with constant 1^ impurity potential has been 
frequently used [27, 29]. Here, we discuss both the use- 
fulness and limitations of the latter model, in which the 
impurity potential in the band-basis (a, b) is given as 

fb = ( h K\ 

then wc obtain 



If we assume Ih and IL are constant, 



t. 



constant — Ik 



|T aa | 2 l + ^ 2 7V(0) 2 / b 2 (l- M 2 ) 



:, (17) 



where T a t, (T aa ) is inter-band (intra-band) T-matrix. 
Here, u = I^/Ib, and iV(0) is the density-of-states for 
each band. In the Born regime ttN(0)I^ <C 1, then 



t. 



constant — 7k 



u holds, and therefore the s±-wave state 



is fragile against impurities except when ti<l. 
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In the unitary regime nN(0)Ih ^ 1, the s±-wave 
state is robust against impurities unless |u| = 1 since 
eq. (17) decreases in proportion to I^ 2 . However, this 
result is totally changed in the five-orbital model, in 
which I h (k, k') = U^(k) ■ W(k, k') ■ U(k') is momentum- 
dependent: Once lb is fc-dependent, then eq. (17) does 
not hold as proved in Ref. [30]. Instead, the relation 

*i b (M0 ~ " 2 ( 18 ) 

holds for all the impurity models studied here even in the 
unitary regime, as recognized by the numerical analysis 
in Fig. 9. Therefore, the constant /j, model is applicable 
to Fe-based superconductors only for Born impurities. 



case would be the lower limit for realistic impurity mod- 
els for Fe-based superconductors. 

In this paper, we studied the effect of in-plane impuri- 
ties. In the case of "out-of-plane" impurities, the radius 
of the impurity potential <f> could be much longer than 
the lattice spacing aFc-Fc [48]. In this case, the relation 
U = \I' b /h\ <C 1 might be realized because of the suppres- 
sion of large angle scattering, and then the s±-wave state 
would be robust against impurities. Therefore, study of 
the out-of-plane impurities would not be useful to distin- 
guish between s ++ - and s±-wave states. 



VI. SUMMARY 



0.2 0.4 0.6 0.8 1 




0.2 0.4 ,0.6 0.8 1 



FIG. 10: (Color online) Obtained -AT c /y of the s±-wave 
state in the two-band constant lb model studied in Ref. [27]. 
The horizontal axis is x = 2u 2 /(l + u 2 ) = 2I h 2 /(I h 2 + I 2 ). 
J b — lb corresponds to x — 1. 

Figure 10 shows the — AT c /7 in the two-band constant 
l b model as function oix = 2u 2 /{l+u 2 ) = II' 2 / \I h 2 +I 2 ), 
in case that the s±-wave state due to inter-band repulsion 
is realized. In this model, —AT c /j is given by eq. (12) of 
Ref. [27]. Note that po oc 7. We see that the relationship 
— AT c /7 oc x holds only in the Born regime nN(0)Ib <S 
1. However, the relation —AT c /j <C x is realized in the 
unitary and intermediate regimes except for |zt| ~ 1. 

Based on the constant lb-model, Ref. [32] mentioned 
that the s±-wave state with T c q = 30K disappears for 
p c T « 100 (1000) [//ficm] at u — 0.5 (0.2) in the inter- 
mediate regime 7riV(0)/b ~ 2.5, and tried to explain the 
experimental small impurity effect on T c based on the 
s_l--wave scenario assuming that u -C 0.5. However, we 
obtain t = u ~ 1 for the realistic impurity models of 
3d-impurity atoms as shown in Fig. 9. In addition, p" 
given in Ref. [32] is very overestimated for u <C 0.5, 
since they apply the constant Jb-model to the interme- 
diate regime and the CVC for the resistivity is dropped. 
In fact, we obtain p" = 17.3z _1 ~ 50 [/jficm] for the 
impurity-induced NL-00 in Fig. 8; t = u ~ 0.4 in this 



Based on the realistic five-orbital model for Fe-based 
superconductors, we have presented a systematic study 
of the impurity effect on the s±-wave SC state. We stud- 
ied impurity model parameters for 3d- and 4d-impurity 
atoms obtained by a recent first principle study [33], in 
addition to various non-local impurity models. The ob- 
tained values of R = — AT c /po for 7ii mp <C 1 as function 
of x = 2|r in tor| 2 /(|Tl n tor| 2 + [TWal 2 ) are summarized in 
Fig. 9. According to the Abrikosov-Gorkov theory, R is 
essentially independent of T c q. 

The main results are shortly summarized as follows: 

(i) For 3e?-impurity atoms, |T intcr |/|r intra | > 0.9 is 
realized, and we obtain R = 3z ~ 5z [K//xQcm], which 
is comparable to Ri =rx > = 3.6z [K//xOcm] for the on-site 
model with I = 00. Also, p" ~ 5z _1 [/ificm] for 
T c0 = 30K. 

(ii) For Ru-impurity atoms, |Ti n ter|/|?intra| = 0.6 — 0.7 
is realized, and both R and 1 / p" are about half of those 
for 3d-impurity atoms. 

(iii) For impurity-induced NL-00 model, which 
gives a very wide-range impurity potential, we obtain 
iTintorl/lTlntral ~ 0.4 and R - 1.5z [K/^cm]. Similar 
result is obtained for the VA s -model and ^i-model with 
finite /. 

(iv) The CVC is important to obtain reliable R 
and p": R is approximately doubled by CVC for Vi- 
and VAs-impurity models, for example. 

(v) In case we use the "constant 1^ impurity model" by 
putting \I h /h\ = |Ji n tcr|/|Tintra|, we have to keep in the 
Born regime nN(0)Ih <C 1, especially when \I b /Ib\ -C 1. 

The abovementioned results are essentially consistent 
with the previous theoretical study for the on-site im- 
purity model by Onari and Kontani [30]. Experimen- 
tally, the critical value of po for the disappearance of 
T c o ~ 30K is 300 ~ 500 [/iSlcm] , which means that the av- 
eraged mean-free-path is comparable to the lattice spac- 
ing (~ 0.3nm) as discussed in Ref. [30]. Therefore, exper- 
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imentally observed robustness of T c against various kinds 
of impurities in optimally-doped Fc-based superconduc- 
tors [22-25] would be consistent with the s++-wave state. 

In this paper, we have introduced only the repul- 
sive pairing interaction due to spin fluctuations. When 
both spin and orbital fluctuations strongly develop, inter- 
orbital attractive interaction and intra-orbital repulsive 
interaction coexist. This situation is naturally expected 
near the orthorhombic phase, and is actually reproduced 
by the SC-VC theory based on the Hubbard model [16]. 
Then, the present study indicates that impurity-induced 
crossover form s±-wave to s ++ -wave states would be re- 



alized in some Fe-based superconductors [13, 49]. 



Acknowledgments 

We are grateful to M. Sato, Y. Kobayashi, J. Li, Y. 
Matsuda, T. Shibauchi, and M. Nakajma for valuable 
discussions. This study has been supported by Grants- 
in-Aid for Scientific Research from MEXT of Japan. Nu- 
merical calculations were partially performed using the 
Yukawa Institute Computer Facility. 



[1] Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono: 

J. Am. Chem. Soc. 130, 3296 (2008). 
[2] K. Hashimoto, et al.: Phys. Rev. Lett. 102 (2009) 

017002. 

[3] M. A. Tanatar, et al.: Phys. Rev. Lett. 104 (2010) 
067002. 

[4] P. J. Hirschfeld, M. M. Korshunov, and I. I. Mazin, Rep. 

Prog. Phys. 74 (2011) 124508. 
[5] M. Yi, D. H. Lu, J.-H. Chu, J. G. Analytis, A. P. Sorini, 

A. F. Kemper, S.-K. Mo, R. G. Moore, M. Hashimoto, 

W. S. Lee, Z. Hussain, T. P. Devereaux, I. R. Fisher, 

Z.-X. Shen, PNAS 108, 6878 (2011). 
[6] R.M. Fernandes, L. H. VanBebber, S. Bhattacharya, P. 

Chandra, V. Keppens, D. Mandrus, M.A. McGuire, B.C. 

Sales, A.S. Sefat, and J. Schmalian, Phys. Rev. Lett. 105, 

157003 (2010) 

[7] M. Yoshizawa et al, Phys. Soc. Jpn. 81, 024604 (2012). 
[8] T. Goto et al, J. Phys. Soc. Jpn. 80, 073702 (2011). 
[9] K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, H. 

Kontani, and H. Aoki, Phys. Rev. Lett. 101, 087004 

(2008). 

[10] I. I. Mazin, D.J. Singh, M.D. Johannes, and M.H. Du: 

Phys. Rev. Lett. 101 (2008) 057003. 
[11] S. Graser, G. R. Boyd, C. Cao, H.-P. Cheng, P. J. 

Hirschfeld, and D. J. Scalapino, Phys. Rev. B 77, 

180514(R) (2008). 
[12] A. V. Chubukov, D. V. Efremov, and I. Eremin, Phys. 

Rev. B 78, 134512 (2008). 
[13] H. Kontani and S. Onari, Phys. Rev. Lett. 104, 157001 

(2010). 

[14] T. Saito, S. Onari, and H. Kontani, Phys. Rev. B 82, 
144510 (2010). 

[15] T. Saito, S. Onari, and H. Kontani, Phys. Rev. B 83, 

140512(R) (2011) 
[16] S. Onari and H. Kontani, Phys. Rev. Lett. 109, 137001 

(2012) 

[17] Y. Ohno, M. Tsuchiizu, S. Onari and H. Kontani, J. 
Phys. Soc. Jpn. 82 (2013) 013707 

[18] T. Hanaguri, S. Niitaka, K. Kuroki, and H. Takagi, Sci- 
ence 328, 474 (2010) 

[19] A. D. Christianson, E. A. Goremychkin, R. Osborn, 
S. Rosenkranz, M. D. Lumsden, C. D. Malliakas, I. S. 
Todorov, H. Claus, D. Y. Chung, M. G. Kanatzidis, R. 
I. Bewley, and T. Guidi, Nature 456, 930 (2008). 

[20] D. S. Inosov, J. T. Park, P. Bourges, D. L. Sun, Y. 
Sidis, A. Schneidewind, K. Hradil, D. Haug, C. T. Lin, 



B. Keimer, and V. Hinkov .Nature Physics 6, 178 (2010) 
[21] S. Tatcmatsu, Y. Yasui, T. Moyoshi, K. Motoya, K. 
Kakurai, and M. Sato, to be published in J. Phys. Soc. 
Jpn. 

[22] A. Kawabata, S. C. Lee, T. Moyoshi, Y. Kobayashi, and 

M. Sato, J. Phys. Soc. Jpn. 77 (2008) Suppl. C 103704; 

M. Sato, Y. Kobayashi, S. C. Lee, H. Takahashi, E. 

Satomi, and Y. Miura, J. Phys. Soc. Jpn. 79 (2009) 

014710; S. C. Lee, E. Satomi, Y. Kobayashi, and M. Sato, 

J. Phys. Soc. Jpn. 79 (2010) 023702. 
[23] Y. Nakajima, T. Taen, Y. Tsuchiya, T. Tamegai, H. Kita- 

mura, and T. Murakami, Phys. Rev. B 82, 220504 (2010). 
[24] J. Li, Y. F. Guo, S. B. Zhang, J. Yuan, Y. Tsujimoto, X. 

Wang, C. I. Sathish, Y. Sun, S. Yu, W. Yi, K. Yamaura, 

E. Takayama-Muromachiu, Y. Shirako, M. Akaogi, and 

H. Kontani, Phys. Rev. B 85, 214509 (2012) 
[25] K. Kirshenbaum, S. R. Saha, S. Ziemak, T. Drye, and J. 

Paglione, Phys. Rev. B 86, 140505 (2012) 
[26] S. Onari, H. Kontani, and M. Sato, Phys. Rev. B 81, 

060504(R) (2010); S. Onari and H. Kontani, Phys. Rev. 

B 84, 144518 (2011) 
[27] Y. Senga and H. Kontani: J. Phys. Soc. Jpn. 77 (2008) 

113710; 

[28] Y. Bang, HY. Choi and H. Won, Phys. Rev. B 79, 
054529 (2009). 

[29] Y. Senga and H. Kontani, New J. Phys. 11, 035005 
(2009). 

[30] S. Onari and H. Kontani, Phys. Rev. Lett. 103 (2009) 
177001. 

[31] Z.J. Yao, W.Q. Chen, Y.K. Li, GH. Cao, H.M. Jiang, 
Q.E. Wang, Z.A. Xu and F.C. Zhang, Phys. Rev. B 86, 
184515 (2012) 

[32] Y. Wang, A. Kreisel, P. J. Hirschfeld and V. Mishra, 

arXiv: 1210.7474 
[33] K. Nakamura, R. Arita and H. Ikeda, Phys. Rev. B 83, 

144512 (2011) 

[34] T. Kawamata, E. Satomi, Y. Kobayashi, M. Itoh, and M. 

Sato, J. Phys. Soc. Jpn. 80 (2011) 084720 
[35] P. C. Canfield, S. L. Budko, Ni Ni, J. Q. Yan, and A. 

Kracher, Phys. Rev. B 80, 060501(R) (2009). 
[36] N. Ni, A. Thaler, J.Q. Yan, A. Kracher, E. Colombier, S. 

L. Budko, P. C. Canfield, and S. T. Hannahs, Phys. Rev. 

B 82, 024519 (2010). 
[37] S. Ideta, T. Yoshida, I. Nishi, A. Fujimori, Y. Kotani, 

K. Ono, Y. Nakashima, S. Yamaichi, T. Sasagawa, M. 

Nakajima, K. Kihou, Y. Tomioka, C. H. Lee, A. Iyo, H. 



11 



Eisaki, T. Ito, S. Uchida and R. Arita, Phys. Rev. Lett. 

110, 107007 (2013) 
[38] H. Kontani, Rep. Prog. Phys. 71, 026501 (2008) 
[39] Y. Inoue, Y. Yamakawa and H. Kontani, Phys. Rev. B 

85, 224506 (2012) 
[40] H. Ikeda, K. Nakamura and R. Arita, unpublished. 
[41] H. Wadati, I. Elfimov and G. A. Sawatzky, Phys. Rev. 

Lett 105, 157004 (2010) 
[42] A. F. Kemper, C. Cao, P. J. Hirschfeld, and H.-P. Cheng, 

Phys. Rev. B 80, 104511 (2009); ibid, Phys. Rev. B 81, 

229902(E). 

[43] J.-H. Chu, J. G. Analytis, K. D. Greve, P. L. McMahon, 
Z. Islam, Y. Yamamoto, and I. R. Fisher, Science 329, 
824 (2010); J. J. Ying, X. F. Wang, T. Wu, Z. J. Xiang, 
R. H. Liu, Y. J. Yan, A. F. Wang, M. Zhang, G. J. Ye, P. 
Cheng, J. P. Hu, and X. H. Chen, Phys. Rev. Lett. 107, 



067001 (2011) 

[44] S. Kasahara, H.J. Shi, K. Hashimoto, S. Tonegawa, Y. 

Mizukami, T. Shibauchi, K. Sugimoto, T. Fukuda, T. 

Terashima, A. H. Nevidomskyy and Y. Matsuda, Nature 

486, 382 (2012) 
[45] T.-M. Chuang, M.P. Allan, J. Lee, Y. Xie, N. Ni, S.L. 

Budko, G.S. Boebinger, P.C. Canfield and J.C. Davis, 

Science 327, 181 (2010) 
[46] C.-L. Song et al, Science 332, 1410 (2011). 
[47] H. Kontani and M. Ohno, Phys. Rev. B 74, 014406 (2006) 
[48] S. Graser,P. J. Hirschfeld, LY. Zhu and T. Dahm, Phys. 

Rev. B 76, 054516 (2007). 
[49] D. V. Efremov, M. M. Korshunov, O. V. Dolgov, A. 

A. Golubov and P. J. Hirschfeld, Phys. Rev. B 84, 

180512(R) (2011) 



